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ABSTRACT 

The ultraluminous compact X-ray sources (ULXs) generally show a curving spec- 
trum in the 0.7-10 keV ASCA bandpass, which looks like a high temperature ana- 
logue of the disk dominated high/soft state spectra seen in Galactic black hole binaries 
(BHBs) at high mass accretion rates. Several ULXs have been seen to vary, and to 
make a transition at their lowest luminosity to a spectrum which looks more like a 
power law. These have been previously interpreted as the analogue of the power law 
dominated low/hard state in Galactic BHBs. However, the ULX luminosity at which 
the transition occurs must be at least 10-50 per cent of the Eddington limit assum- 
ing that their highest luminosity phase corresponds to the Eddington limit, while for 
the Galactic BHBs the high/soft-low/hard transition occurs at a few per cent of the 
Eddington limit. Here we show that the apparently power law spectrum in a ULX in 
IC342 can be equally well fit over the ASCA bandpass by a strongly Comptonised opti- 
cally thick accretion disk with the maximum temperature of ~ 1 keV. Recent work on 
the Galactic BHBs has increasingly shown that such components are common at high 
mass accretion rates, and that this often characterises the very high (or anomalous) 
state. Thus we propose that the power law type ULX spectra are not to be identified 
with the low/hard state, but rather represent the Comptonisation dominated very 
high/ anomalous state in the Galactic BHBs. 

Key words: galaxies: spiral - galaxies: individual: IC 342 - X-rays: galaxies - X-rays: 
binaries 



1 INTRODUCTION 

Since the Einstein era, many ultraluminous compact X-ray 
sources (ULXs; Makishima et al. 2000) with X-ray luminosi- 
ties Lx > 10 39-40 erg s _1 have been found in spiral arms of 
nearby spiral galaxies (e.g., Fabbiano 1989). Such high val- 
ues of Lx exceed the Eddington limit, Le, for a neutron star 
by several orders of magnitude, and instead suggest massive 
(30-100 Mq) accreting black holes (BHs). However, it is 
difficult to understand how such massive BHs could form: 
single massive stars have extreme mass loss throughout their 
life, and the maximum BH mass expected is of order 10-15 
Mq (e.g. Fryer & Kalogera 2001). Since there are no nearby 
systems either in our galaxy or in M31 which could be easily 
studied then the nature of the ULX remained a mystery. 

ASCA (Tanaka, Inoue, & Holt 1994) data led to a break- 
through for these objects by providing the first moderate res- 
olution energy spectra. As reported by many authors (e.g., 
Makishima et al. 2000, and references therein), it is clear 



that the majority of the most luminous ULXs exhibit spec- 
tra which are well fitted by the multicolor disk model (MCD 
model; Mitsuda et al. 1984), similar to the case of Galac- 
tic/Magellanic BH binaries (BHBs) at high mass accretion 
rates, M (e.g. Makishima et al. 1986; the review by Tanaka 
& Lewin 1995). This, together with the variability of these 
systems including periodic variation (Bauer et al. 2000; Sug- 
iho et al. 2001) and their general association with regions of 
ongoing star formation (Zezas, Georgantopoulos, & Ward 
1999; Roberts & Warwick 2000; Fabbiano, Zezas, & Murray 
2001) has led to their identification with BHBs. 

There are then three possibilities, first that these really 
are intermediate mass BHBs, formed perhaps via mergers of 
massive stars/BHs in a compact star cluster (e.g. Ebisuzaki 
et al. 2001). Alternatively, they could be "normal" mass 
(~ IOMq) stellar BHs accreting beyond the critical accre- 
tion rate at which the disk luminosity Ldisk reaches Le- Re- 
cent work suggests that the Eddington limit can be violated 
by the disk becoming clumpy (e.g., Krolik 1998; Gammie 
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1998; Begelman 2002 and references therein), and it has 
long been known observationally that super critical accre- 
tion can happen (e.g. Cir X-l). A third alternative is that 
these are "normal" mass BHs, but that their X-ray luminos- 
ity is strongly anisotropic (beamed), so that the bolometric 
luminosity is overestimated. This beaming is highly unlikely 
to be the relativistic beaming seen in jet sources as this 
generally leads to a power-law (PL) type of spectrum (e.g. 
blazars), very unlike the curving spectrum seen by ASCA for 
the majority of the ULX. However, strong anisotropy of the 
disk flux might be produced if the disk is geometrically thick 
in its inner regions. The radiation can then be strongly col- 
limated by a funnel (e.g. Madau 1988), although the factor 
of 10-100 required (King et al. 2001) seems extreme. Such 
an extreme thick disk might be expected to form only under 
a super critical accretion rate. So this is merely an adden- 
dum to the high M scenario, rather than an independent 
alternative. 

All these alternatives involve an extreme of one kind 
- either mass, radiation luminosity or disk shape. One way 
to test these is to compare the ULX spectra and spectral 
variability with those of Galactic BHBs. This has become 
much more feasible in recent years with the unprecedented 
volume of data from the Galactic BHBs gathered by RXTE. 
Here we use the bright Galactic BHB transient XTE J1550 — 
564 to observationally determine the spectra and spectral 
variability of high M disks, and show that the ULX are 
indeed compatible with being massive (30-100 Mq) BHs 
accreting at close to the critical accretion rate. 



2 A COMPARISON OF THE ULX WITH BHB 

Here we specify the puzzles in understanding ULXs as 
~ 100 Mq BHs. It has long been recognized that the Galac- 
tic/Magellanic BHBs reside in either of the two distinct spec- 
tral states, the high/soft state or the low/hard state. In the 
hard state, in which M is generally low, the BHBs exhibit 
a hard PL spectrum in the 2-20 keV band. By contrast, 
the soft state is generally seen when M is high, and it is 
characterized by the MCD spectrum, which approximates 
the optically thick standard accretion disk (Shakura & Sun- 
yaev 1973). Based on this viewpoint, the MCD-type ULXs 
have so far been regarded as residing in the usual soft state, 
even putting aside the issue of how to make the required 
~ 100 M Q BHs. 

An increasing number of ULXs shows spectral transi- 
tions from an MCD (curving) type of spectrum to one which 
is better described by a hard PL. This has been interpreted 
as the transition to the usual hard state (e.g. Kubota et al. 
2000, hereafter Paper I; La Parola et al, 2001; Mizuno 2000). 
However, such a straightforward analogy between ULXs and 
BHBs gives the following puzzles. 

(i) The "soft state" ULXs have inner disk temperatures 
Tin which are too high for their implied high masses. Equiv- 
alently, when the BH mass is estimated from the Eddington 
argument, R 1Yl with reasonable correction for a boundary 
condition (Kubota et al. 1998) and spectral hardening factor 
(Shimura & Takahara 1995) falls much below the last stable 
orbit for a non-spinning BH, 6i? g , where R g — GM/c 2 is the 
gravitational radius. This contrasts with soft-state BHBs, 
where Ri n generally agrees with 6R g . 



(ii) The value of Ri n is time variable, in five bright "soft 
state" ULXs including the two in IC 342 (Mizuno et al. 
2001). This again makes a contrast to the case of soft-state 
BHBs, where R ln is approximately constant for each source. 

(iii) Assuming that the "soft state" corresponds to Le, 
the threshold luminosity for the spectral transition of ULXs 
(~ 0.3 Le; Paper I) becomes much higher than is seen among 
BHBs, typically 0.01-0.03 L E . Of course, the MCD-statc 
luminosity could be much below Le , but then the BH mass 
in ULXs would have to be even higher. 

Makishima et al. (2000) attribute the first puzzle to an 
extreme BH rotation, because the last stable orbit is then 
reduced to 1.23i? g , although the other two problems remain 
unsolved. They suggested the scenario of Kerr BHs shining 
at Le for ULXs. Following Makishima et al. (2000), Ebisawa 
et al. (2001) showed that, with standard Shakura-Sunyaev 
disks, even with Kerr BHs it is difficult to obtain the required 
high temperatures without super (or near) critical accretion 
rates. A more plausible explanation is that at high M the 
disk structure changes from that of the standard disk due 
to the disk being so optically thick that radial advection of 
the radiation becomes important. These disks, called "slim 
disk" (Abramowicz et al. 1988), are different in structure 
from the standard disks. Pressure support becomes impor- 
tant, so the material is not in Keplarian rotation. The inner 
edge of the optically thick accretion disk is then not neces- 
sarily at the last stable orbit, but can be closer to the BH 
(e.g. Abramowicz et al. 1988; Watarai et al. 2000). The small 
and changing inner disk radius (puzzles i and ii) could be 
explained in a Schwarzschild metric if the disk penetrates in- 
creasingly further into the plunging region at high M rather 
than abruptly truncating at the last stable orbit, so Mizuno 
et al. (2001) and Watarai, Mizuno & Mineshige (2001) pro- 
posed that MCD-type ULXs may have slim accretion disks 
rather than standard ones. 

The inner disk radius, however, is very dependent on the 
viscosity prescription. For small viscosity it can decrease to 
~ 47? g for super Eddington luminosities, whereas for large 
viscosity it can be slightly larger than 6i? g (e.g. Abramow- 
icz et al. 1988; Artemova et al., 2001). Also, this mecha- 
nism for producing a decreasing radii can only be used for 
Schwarzschild BHs as in extreme Kerr BHs there is so lit- 
tle space between the last stable orbit and the horizon. The 
spectra of these slim disks have been calculated by Watarai 
& Mineshige (2001) (Schwarzschild, low viscosity) and Be- 
loborodov (1998) (Schwarzschild and Kerr, high viscosity). 
Since some of the emission from the smallest radii is ad- 
verted rather than radiated, the spectra have less high tem- 
perature emission than a standard disk, but this can be 
somewhat compensated by the decrease in inner radius for 
the Schwarzschild, low viscosity disks. Since the standard 
disk structure calculations have difficultly in producing the 
high temperatures observed from the ULX even with Kerr 
BHs for sub-critical rates (Ebisawa et al. 2001), then it seems 
unlikely that the slim disks will substantially help problems 
(i) and (ii) unless the accretion rate is super- critical. 

Therefore it seems likely that the standard curving 
MCD-type ULX spectra are from super critical accretion. In 
which case the third puzzle becomes acute. Recently, Kub- 
ota et al. (2001a,b) has observationally suggested a novel un- 
derstanding of high-luminosity accretion disks, using RXTE 
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data on some Galactic BHBs. They have identified three dis- 
tinct spectral regimes of the optically-thick accretion disk in 
XTE J1550-564 and GRO J1655-40. Figure 1 shows their 
result on XTE J 1550 — 564 (the distance and inclination an- 
gle are assumed to be 5 kpc and 60°, respectively). 

(1) The standard regime where Ri n remains constant 
when fit with MCD models, i.e. Ldisk oc Tj„. 

(2) The luminous apparently standard regime where the 
disk luminosity rises slightly less quickly with temperature, 
Ldisk oc Tjn 2 , and the spectral shape is slightly distorted 
from the standard-disk one. In the literature, both this and 
(1) are identified as the soft state, with disk dominated spec- 
tra. 

(3) The intermediate anomalous regime, where the disk 
Comptonisation suddenly increases and the spectrum be- 
comes much harder. This is often termed the very high state 
in the literature (Miyamoto et al. 1991). 

Kubota et al. (2001a,b) suggested a picture for this, in which 
the standard regime is where the accretion flow is described 
by the Shakura-Sunyaev geometrically thin disk approxi- 
mation. As the disk luminosity increases, this approxima- 
tion breaks down, and the disk becomes slim. This theo- 
retical picture explains both the change in spectral shape 
(as the hotter, inner disk emission is preferentially advected 
rather than radiated) and the different Ti n -Ldiak relationship 
(Watarai et al. 2000) seen in the apparently standard regime, 
although for the likely mass/distance of XTE J1550 — 564 
then the transition threshold appears at O.ILe, rather than 
at the predicted Le- It should be noted that, though Comp- 
tonisation in the disk can shift the apparent disk tempera- 
ture/radius as measured by MCD models, the effect of this 
is Jess important at high M than at low M (Merloni, Fabian 
& Ross 2000), opposite to the observed behaviour on the 
Tln-^disk plane. 

The clear observational result is that the high M Galac- 
tic BHBs show a variety of spectral shapes, which corre- 
late fairly well with L^^/Le- Overlaid on the BHB data in 
Fig. |l| are the results from fitting several spectra from the 
ULX IC 342 source 1 with the MCD model (as in Mizuno 
et al. 2001, except that an inclination angle of 60° is as- 
sumed) . We thus notice a clear similarity between the ULX 
behaviour and that of XTE J 1550 — 564 in its apparently 
standard regime. This similarity argues for a slim disk inter- 
pretation of ULX spectra, as first suggested by Mizuno et 
al. (2001). As to the MCD-type spectrum of ULXs, we here 
simply point out its similarity to the apparently standard 
regime of BHBs, although to get the same temperature from 
a much more massive BH requires that the disk should be 
accreting at a correspondingly larger fraction of the critical 
rate. Overlaid on Fig.l is the standard disk Ti n -Ldisk rela- 
tion for 10, 30 and 100 solar mass Schwarzchild and Kerr 
(a=0.55) black holes. 

In Fig.l, the horizontal solid line indicates the 0.7-10 
keV luminosity of IC 342 source 1 in 2000, when it made 
a spectral transition into the PL-type spectrum. The ex- 
trapolated 1-30 keV luminosity in the PL-type spectrum is 
close to Ldisk in the dimmest MCD spectral state (the 1993 
data). This illustrates the problem (iii) above; if we require 
that the ULX is shining at close to Le, then this PL-type 
state also has a luminosity close to Le, so it is unlike a usual 
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Figure 1. The calculated Ldisk plotted against the observed T ln . 
As for XTE J1550-564 (taken from Kubota 2001a), three kinds of 
symbols specify the data obtained in the standard regime (open 
circles), the apparently standard regime (filled circles), and the 
anomalous regime (gray crosses). For the anomalous regime data, 
^disk + Lthc are plotted instead of Ldisk- The results of IC 342 
source 1 in 1993, the MCD-type state, are also plotted with open 
squares. For comparison, the 0.7-10 keV absorption corrected lu- 
minosity of the same source in the PL-type state is presented as a 
horizontal solid line. The bolometric luminosity calculated under 
the thcomp model is also plotted against T; n with a filled gray 
square. The dash-dotted lines represent the standard Ldisk oc TJ* 
relation for Schwarzshild BHs of 10, 30, and 100 Mq, while the 
dashed lines represent the same relation for mildly rotating BHs 
(a = 0.55). Relativistic effects are not taken into account to calcu- 
late these lines. The dotted lines represent Ldisk °c T? relation. 
The Eddington limits are also shown for Shwarzshild BHs and 
mild Kerr BHs with a dash-dotted line and a dashed line, respec- 
tively. 

low/hard state for the stellar BHBs, even though it has a 
characteristic PL spectral shape of the low/hard state. 

Here we propose instead that the PL-type spectrum 
seen in the ULX marks a transition to the Comptonised 
anomalous regime rather than the usual low/hard state. 



3 REANALYSIS OF PL- TYPE SPECTRUM 

In these few years, our understanding of both ULXs and 
stellar BHBs has thus made rapid progress observationally 
and theoretically. This condition makes us possible to re- 
consider the PL-type ULXs and solve the third puzzle pre- 
sented in §2. Accordingly, we re-examine the GIS data of 
IC 342 source 1 in 2000, the longest observation among the 
PL-type ULXs. The observational details are given in Paper 
I and Sugiho et al. (2001). 

3.1 Characteristics of the spectrum 

The overall GIS spectrum is well described by an absorbed 
PL of r = 1.73 ± 0.06 (Paper I). In order to investigate 
whether this is truly a simple PL shape, we fit the low and 
high energy band separately. The 0.7-4 keV GIS data give 
a flatter power law index, V — 1.54 ± 0.12, modified by a 
low energy absorption of Nn = 5.2 ± 0.3 x 10 21 cm -1 , than 
that obtained from fitting the 5-10 keV data with a single 
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Figure 2. The GIS spectrum of IC 342 source 1 in 2000, the 
PL-type state, with a prediction of the best-fit model obtained 
by 0.7-4 keV spectral fitting. A confidence contour of V against 
Nn is also shown. The bottom panel indicate the fit residuals. 

PL where P = 2.4 ± 0.3. Figure ^ shows the extrapolation 
of the 0.7-4 keV spectral fit to the full GIS energy band, 
clearly showing that the spectrum is slightly curved. This 
could be due to spectral complexity around the iron line 
and edge energies as suggested in Paper I. However, the 
low/hard state PL-spectra of galactic BHB usually show an 
opposite curvature, with the 5-10 keV spectrum going harder 
than the 0.7-4 keV spectrum due to frequently observed 
soft excess softening the low energy spectrum, and some 
reflection hardening the higher energy bandpass. 

Thus the subtle downwards curvature of the spectrum 
makes it look different from the BHB low/hard state. In 
order to characterize the spectral shape from another view- 
point, we compare the PL-type spectrum in 2000 with the 
faintest MCD-type spectra seen from the same source in 
1993 (using the same data selection as Mizuno et al. 2001). 
As shown in table 1, the faintest 1993 GIS data is well fit 
with the MCD model of T in = 1.3 keV modified by low 
energy absorption column of A^h = 3.4 x 10 21 cm -2 . These 
parameters are consistent with those in Mizuno et al. (2001). 
Figure |^ shows the ratio of the PL data in 2000 with the best 
fit MCD model to the 1993 faintest data. This is remarkably 
flat below 3 keV, showing that the shape of the low energy 
spectrum does not change. A better explanation of the PL- 
type spectrum could be that it is not a simple PL but has 
a soft component which is a fainter version of that in the 
faintest MCD state, together with a Comptonised tail, sim- 
ilar to the anomalous type (or very high state) spectra seen 
in XTE J1550 - 564 and GRO J1655 - 40. 

3.2 Fit with the Comptonised model 

We re-fit the same PL-type GIS spectrum incorporating 
Comptonisation. We use the same model as for XTE J1550 — 
564 in the anomalous regime (Kubota 2001a), i.e. an approx- 
imate model for thermal Comptonisation based on the Kom- 
paneets equation (thcomp: Zdziarski Johnson & Magdziarz 
1996, Zycki, Done, & Smith 1999). We choose a disk black- 
body as the seed photon distribution, and include a sep- 
arate uncomptonised MCD whose temperature is tied to 
that of the seed photons, as the thcomp model only cal- 
culates Comptonised component. We do not include reflec- 
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Figure 3. The PL-type spectrum of IC 342 source 1 (in 2000) 
normalized to the best-fit spectrum of the same source in the 
MCD-type state (in 1993). The latter is convolved with the in- 
strumental response calculated for the position of the source in 
2000. 

tion as the statistics are not sufficient to constrain the fit. 
The fit parameters are maximum temperature of the disk 
blackbody seed photon spectrum, Tin, the electron temper- 
ature T e , photon index T th which expresses spectral shape of 
the thcomp below T e , and the normalizations of the Comp- 
tonised and uncomptonised disk photons. We could not con- 
strain the values of T c with the 0.7-10 keV GIS data. There- 
fore, we fix T c at 20 keV as seen in XTE J 1550 -564 (Wilson 
& Done 2001). 

We show the best fit thcomp parameters in table 1. The 
thcomp model with a low energy absorption can successfully 
reproduce the observed spectrum, with xVdof = 102/86. 
The obtained value of Tin = 1.1 keV is slightly below the 
temperature of 1.3 keV seen in the dimmest MCD-type spec- 
trum in 1993. Although the direct MCD component was 
added to the spectral model, the data do not require it with 
a 90 % upper limit to the direct disk emission being less 
than 25 % of the total 0.7-10 keV energy flux. The lack of a 
direct MCD component is consistent with the optical depth 
of r ~ 3 inferred from the T th and (fixed) T e . 



4 DISCUSSION 

In §3, we have shown that the apparently PL-type spectra 
seen from IC342 in 2000 is most probably not related to the 
low/hard spectra in Galactic BHBs. It shows significant de- 
viations from a PL shape, in the sense that the spectrum 
softens at higher energies. This is opposite to the behaviour 
of the Galactic BHBs in the low/hard state. We propose that 
the PL-type spectra are instead the analogue of the anoma- 
lous regime (also termed the very high state) spectra seen in 
the Galactic BHs at high luminosities. We demonstrate this 
by fitting Comptonisation models, thcomp, to the PL-type 
data, and show that they can indeed give as good a fit to the 
data as a PL continuum plus ionized edge (Paper I). If the 
ionized edge is not incorporated, the PL fit is significantly 
inferior to the thcomp fit (table 1), because of the intrinsic 
concaveness of the spectral shape. 

Additionally, validity of this interpretation can be 
tested by investigating the location of the obtained result on 
the Ldisk-Tin plane (Fig.l). There, we also show the 0.1-100 
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Table 1. Spectral parameters of IC 342 source 1 with 90% confidence limits 



epoch 


model 


range(keV) 




r or r th 


Tin (keV) 


flux 6 


X 2 /dof 


2000 


PL 
PL 


0.7-4 
5-10 


5.2 ±0.3 


1.54 ±0.12 
2.4 ±0.3 




4.0 
6.1 


71.7/57 
14.4/19 




thcomp 
PL C 
PL d 

MCD e 


0.7-10 
0.7-10 
0.7-10 
0.7-10 


3.2 ±0.4 
6.4 ±0.7 

6.8 ±0.7 

1.9 ±0.4 


2.2 ±0.4 
1.73 ±0.06 
1.78 ±0.06 


1.1 ±0.3 
2.06 ±0.08 


3.3 
3.7 
3.9 
2.1 


102.2/86 
101.1/86 
119.4/88 
120.5/88 


1993 (faintest phase) 


MCD 


0.7-10 


3 4+2.1 




1 30+ 019 
1 - JU -o.i6 


6.4 


23.1/33 



., , .0 21 

6 The 0.7—10 keV source flux at the top of the atmosphere, in units of 10 — 12 erg s — 1 cm -2 

c Taken from Paper I. An ionized Fe-K edge at 8.4 ± 0.3 keV is applied, with an optical depth of 0.9 ± 0.5. 

Fit result by a single PL with absorption column. 
e Taken from Paper I 



keV thcomp luminosity, L tn c = 1.1 X 10 40 erg s" 1 cm" 2 , ex- 
trapolated from the best fit model assuming isotropic emis- 
sion. If the Comptonisation is from overheating of the disk 
(Beloborodov 1998) then this represents the disk luminosity. 
If instead it is from a separate corona, then the disk lumi- 
nosity is amplified by Comptonisation. However, at this low 
compton y parameter the amplification is rather low (less 
than a factor of 2), so we use the thcomp luminosity as our 
estimate for the disk luminosity and plot this with the seed 
photon temperature on Fig [lj. The PL-type spectrum then 
lies nicely on the same luminosity-temperature relation as 
defined by the MCD-type spectra. This predicts a ULX mass 
between 30 M Q (a = 0) and 150 M (a = 0.998) if the PL- 
type spectrum marks the break between the standard and 
apparently- standard (slim disk) regimes. This break then ap- 
pears at ~ 1-2Le- If instead the break is at ~ 0.1-0.21/e as 
in J1550 then this implies a mass of 50-250 Mq. 

Based on these observational results, including spectral 
softening at higher energies, X-ray luminosity, and disk tem- 
perature within the framework of the disk Comptonisation, 
we conclude that the PL-type spectra of IC 342 source 1 in 
2000 are related to the standard MCD-type disk with strong 
disk Comptonisation state rather than to the low/hard state, 
and that the transition between 1993 and 2000 is likely an 
apparently standard to anomalous transition, rather than the 
canonical high/soft to low/hard transition. 
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